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Abstract 

Experimental result regarding the maximum limit of the radius of the electron ~ 10~ 16 cm and a 

few of the theoretical works suggest that there might have some negative energy density regions 

within the particle in general theory of relativity. It is argued in the present investigation that 

such a negative energy density also can be obtained with a better physical interpretation in the 
(N ■ 

■ framework of Einstein-Cartan theory. 

o 

<N ' 1. INTRODUCTION 

t— i ; 

CN , Recently, Cooperstock and Rosen [1989], Bonnor and Cooperstock [1989], and Herrera and 
O ■ 

Varela [1994] have shown that within the experimentally obtained upper limit of the size of 
the electron ( ~ 10 -16 cm ) [Quigg, 1983], when it is modeled as a charged sphere obeying 
(3JT), Einstein-Maxwell theory, must contain some negative gravitational mass density regions within 
the particle . According to Cooperstock, Rosen and Bonnor (CRB) [1989], the rest mass or 
^ ■ active gravitational mass within this sphere, by virtue of the relation 

Q 2 

M = m-—, (1) 

is negative and about 10~ 52 cm (when the inertial mass or effective gravitational mass, charge 
and radius, respectively, of the electron, are m = 6.76 x 10~ 56 cm, q = 1.38 x 10 -34 cm, and 
a = 10~ 16 cm in relativistic units ). Further, Herrera and Varela (HV) [1994] have shown, in 
one of the cases of their paper, that the matter-energy density 

p = (ag 2 + ^ 2 )(a 2 -r 2 ), (2) 

for the constant a = —4.77 x 10 95 cm -6 (when radius a ~ 10~ 16 cm) is also negative, do being 
the constant charge density at the centre of the spherical distribution. These models, however, 
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lack spin and magnetic moment and hence do not possess the actual physical characteristics 
required for an electron. 

As an alternative way both the groups suggest the stationary Kerr-Newman (KN) metric [New- 
man et al., 1965] related to the solution of Einstein-Maxwell equations [Misner et al., 1973] to 
be more appropriate than those described earlier. However, in this context it is also to be men- 
tioned here that the KN metric cannot be valid for distance scales of the radius of a subatomic 
particle [Mann and Morris, 1993; Herrera and Varela, 1994]. 

We, therefore, feel that the problem can be tackled in the framework of Einstein-Cartan (EC) 
theory, where torsion and spin are inherently present in the formulation of the theory itself. 

2. AN OVERVIEW: THE NEGATIVE DENSITY MODELS 

Before going into the Einstein-Cartan theory let us have a birds' eye view of the negative matter- 
energy density models already we have mentioned in the introduction, 
(i) THE COOPERSTOCK-ROSEN-BONNOR (CRB) MODEL 

Cooperstock and Rosen [1989] and Bonnor and Cooperstock [1989] in their papers have shown 
that any spherically symmetric distribution of charged fluid, irrespective of its equation of state, 
whose total mass, radius and charge correspond to the observed values of the electron, must have 
a negative energy distribution (at least for some values of the radial coordinate). Considering a 
static spherically symmetric charge distribution with the line-element 

ds 2 = e 2 ^dt 2 - e 2X ^dr 2 - r 2 (d9 2 + ainHdtf) (3) 

they have argued that when the Einstein-Maxwell equation 

R%-±5\R = 8n(T\ im) +T\ {em) ) (4) 

is written in the form 

e -2A = e 2u = l _]_ r (g + e -(v+V E 2y dr (5) 

r Jo 

and hence is equated with the Reissner-Nordstrom exterior metric on the boundary r = a, which 
as usual gives 

1 _ ^ + £_ = i _ I / a ( 87r/9 + e -^E 2 )r 2 dr. (6) 
a cr a Jo 

Then for the previous specifications of mass, charge and radius of the electron it can be shown 
that 

4 - ^ ~ 2 x 1(T 36 > 0. (7) 
a z a 
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So, the left hand side of the above equation (6) must be greater than unity and hence on the 
right hand side p < for some values of r implying that the electron must contain some negative 
rest mass density though the net mass is as usual a positive quantity 
(ii) THE HERRERA-VARELA (HV) MODEL 

Following the CRB model [1989] Herrera and Varela [1994] have discussed the fact that the 
electron, when modeled as a relativistic spherically symmetric charged distribution of matter, 
must contain some negative rest mass if its radius is not larger than ~ 1CP 16 cm. In this regard 
they have analyzed some extended electron models and have shown that negative energy density 
distributions result from the requirement that the total mass of these models remains constant 
in the limit of a point particle. Among all these extended electron models the model of Tiwari et 
al. [1984] demands special attention to us which will be seen very much relevant to our present 
work. Herrera and Varela [1994] generalize this model of Tiwari et al. [1984] by introducing a 
condition of anisotropy in the form 

PL ~ Pr = aq 2 r 2 (8) 

where a is a constant. 

Thus the solution obtained by Herrera and Varela [1994] is as follows: 

e -2A = e 2„ = 1 _ H^V^a 2 - 2r 2 ) - -naq 2 r 2 (ba 2 - 3r 2 ), (9) 
45 15 

Pr = -(aq 2 + ^7ra 2 ){a 2 -r 2 ), (10) 

p± = aq 2 r 2 - (aq 2 + ^a 2 ){a 2 - r 2 ), (11) 

m = ^ir 2 a 2 a 5 + ^-iraq 2 a 5 , (12) 

p=(aq 2 + ^a 2 )(a 2 -r 2 ). (13) 

q = ^7TO- a 3 (14) 

The value of a can be obtained from the equation (12) as a = —4.77 x 10 95 cm -6 and hence 
the energy density, as given by the equation (13) is negative for the radius of the electron 
a = 10 -16 cm. 

Now, from the equation (12) it can be seen that the effective gravitational mass, m, is of purely 
electromagnetic origin and corresponds to the TRK model [1984] with a = case. This type 
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of models where mass, including all the other physical parameters, originates from the electro- 
magnetic field alone are known as the electromagnetic mass models [EMMM] in the literature 
[Feynman et al., 1964] and have been investigated by several authors [Florides, 1962, 1983; 
Cooperstock and de la Cruz, 1978; Tiwari et al., 1984, 1986, 1991, 2000; Gautreau, 1985; Gr0n, 
1985, 1986a, 1986b; de Leon, 1987a, 1987b, 1988; Tiwari and Ray, 1991a, 1991b, 1997; Ray et 
al., 1993; Ray and Ray, 1993]. In the present paper we shall construct such a model within the 
framework of Einstein-Cartan theory with negative matter-energy density for some values of the 
radial coordinate. 

3. THE FIELD EQUATIONS OF EINSTEIN-CARTAN THEORY 
The EC field equations are given by 

K i -\PjR = -Kt i i, (15) 

Q'ik - VjQ l ik - $\Q l ji = -KS* jk , (16) 

where t l j is the canonical energy-momentum tensor (asymmetric) , Q % j k is the torsion tensor and 
S l jk is the spin tensor (with k = —8tt, G and c being unity in relativistic units ). 
The asymmetric energy-momentum tensor here is given by 

* i i = T i i + ^V B (S m i *) 1 (17) 

V m being covariant derivative with respect to the torsionless, symmetric Levi-Civita connection 
r l jfc and the symmetric energy-momentum tensor T l j will consist of two parts, viz., matter and 
electromagnetic tensors and which, respectively, are 

T/ m) = (p + p)u\-pg* J , (18) 

r , >m) = ±_ { _ FjkF ik + ±5^F kl F kl ), (19) 

where p is the matter-energy density, p is the fluid pressure, u l is the velocity four-vector (with 
u l Ui = 1) and is the electromagnetic field tensor. 

The conservation equations for the EC theory can be given through the Bianchi identities as 

V fc [(p + P )u k - g ki u l V m (u m S h )] = vPVjp, (20) 
[{p + p)u k - g ki u l V m (u m S li )}V kUj = -Vi(u l Uj ) + u k S jm R m k - \u k S lm R lm 3k (21) 
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Now, electromagnetic fields not being coupled with torsion [Novello, 1976; Raychaudhuri, 1979] 
the Maxwell equations as usual take the form 



VjF*' = J\ (22) 

(^V=^),i = 0. (23) 

The electromagnetic field tensor, Fij, in the above equation (22) is related to the electromagnetic 
potentials as = Aij — Ajj which is equivalent to F^j^ = 0, Ai being the electrostatic 
potentials. Here and in what follows a comma denotes the partial derivative with respect to the 
coordinate indices involving the index. 

Again, the spin tensor and the intrinsic angular momentum density tensor are related as 

S l jk = u l S jk , (24) 

with 

u'Sik = 0. (25) 

Now, assuming that the spins of the individual charged particles composing the fluid distribution 
are all aligned in the radial directions [Prasanna, 1975; Raychaudhuri, 1979; Tiwari and Ray, 
1997] and the matter is at rest with respect to the observer, the non-vanishing components of 
the spin tensor can be obtained, from equations (24) and (25), as 

S°2 3 = S° 32 = S (5oo)" 1/2 , (26) 
whereas, from equation (16), we have the torsion tensor as 

Q°23 = -Q%2 = -««(9do)- 1/2 , (27) 

s = S23 being the only non vanishing component of the intrinsic angular momentum density 
tensor. Here, we have followed the convention (t, r, 8, (f>) = (0, 1, 2, 3). 

The Einstein-Cartan-Maxwell equations with source then can be written as [Tiwari and Ray, 
1997] 

e-^-^) + ^=^~ P + E\ (28) 
e- 2A (^ + l)-l = 8^~ r -£ 2 , (29) 



[z/' + z/ -z/A' + ^ '-]=torp± + E , (30) 
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(r 2 E)' = 47rrVe\ (31) 

where p, p r , p~± and E are the effective matter-energy density, effective pressures (radial 
and tangential) and electric field respectively, and are defined as 

p = p-2irs 2 , (32) 

Pr = Pr — 27TS 2 , (33) 

P±=P±- 2tts 2 , (34) 

E = -exp[-(v + X)]</>' = \, (35) 

Pi Pr, PL, s, <f> an d q being the ordinary matter-energy density, ordinary pressures 
(radial and tangential), spin density, electrostatic potential and electric charge respectively. 
Here, a denotes the electric charge density and prime is used for the derivative with respect to 
the radial coordinate r. 

Also, the conservation equations (20) and (21) here reduce to 

— = -( P + Pr)v + ) + • (36) 

This is the key equation which is to be solved for constructing EMM. 
4. THE SOLUTIONS 

Addition of (28) and (29), under the assumption <7oo<7n = — 1 ( or equivalently, in terms of 
energy-momentum tensors T°o = T l \), provides the pure charge condition 

P + Pr = 0, (37) 

where, in general, p is assumed to be positive and hence p r is negative. However, as is evident 
from equation (32), p, being the effective energy-density, can even be negative due to the positive 
second term related to the spin on the right hand side. Thus, the possibility of equation (33) 
being satisfied with p r being positive is not ruled out. 
To make (31) and (36) solvable, we further assume that 

ae x = a , (38) 

P± ~ Pr = P± ~ Pr = aq 2 r 2 (39) 
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following TRK Model [Tiwari et al., 1984] and HV model [Herrera and Varela, 1994] respectively, 
where do and a are two constants as mentioned earlier. 
By substituting (37) - (39) in (31) and (36), we get 

4 

V 



E=\ = -ira r, (40) 



p r = 2ns 2 - (aq 2 + ^ira 2 )(a 2 - r 2 ), (41) 
p ± = 2-ks 2 - aq 2 (a 2 - 2r 2 ) - ^ira 2 (a 2 - r 2 ), (42) 
P = 2tts 2 + (aq 2 + ?-ira 2 )(a 2 -r 2 ). (43) 



The active gravitational mass 



E 2 

M(r) = 4vr / (p - 2ns 2 + —)r 2 dr (44) 

87T 



takes the form, by virtue of (40) and (43), as 



M(r) = 7rV 2 r 3 [87rcra 6 (5a 2 - 3r 2 ) + 3(5a 2 - 2r 2 )]. (45) 

135 



Thus, the metric potentials A and v are given by 

e _ 2A = e2 , = 1 _2MW ; (46) 



r 



whereas, the effective gravitational mass in (1), can be obtained as 

m = ^7r 2 CT 2 a 5 (l + ^vrcra 6 ), (47) 

which corresponds to the second case (B) of HV model [1994] and of purely electromagnetic 
origin. This again corresponds to the TRK model [1984] with a = case. It can be noted here 
that unlike the matter-energy density the effective gravitational mass is independent of spin. 
In this context it is to be mentioned here that the junction conditions in the EC theory are 
different from that of general theory of relativity and indeed read like this [Arkuszewski et al., 
1975] 

n^|_=0 (48) 

V |-= 2irG(n i S i ) |_ (49) 

where S l is the spin density pseudo-vector. Here condition (48) is the same as in classical 
relativistic hydrodynamics and has already been incorporated by matching the interior solution 



7 



to the exterior Reissner-Nordstrom field at the boundary of the spherical distribution. The 
condition (49), however, is the additional condition to be satisfied in EC theory. In the present 
case, it is only the effective pressure (radial) that vanishes on the boundary and not the ordinary 
radial pressure which, by virtue of equation (41), equals 2tts 2 . The spin is aligned in the 
radial direction and hence the spin density pseudo- vector is hypersurface orthogonal. Thus, the 
boundary condition (49) will become 

Pr\ r =a-0 = 2 ^ s \=a-^ ( 50 ) 

which, depending on whether s is a constant or function of coordinates, will automatically be 
satisfied. 

In this connection it is to be mentioned here that the spin density, s, in the final solutions (41) - 
(43) remains arbitrary (function of r). An explicit functional form of this spin density can also 
be obtained by assuming some additional physically viable possibility, such as the one used by 
Prasanna [1975] by splitting the conservation equation into two parts, the second part relating 
to conservation of spin only, giving the functional form of spin density as s = s§e~ v (where So 
is the value of s at r = i.e. the central spin density). This can, using equations (38) and 
(46) and the condition g 00 gu = 1, (that is, v + A = 0) equivalently be written as so = sqctq = 
constant. The functional form of spin density is, however, not relevant in our discussion as our 
problem is concerned with the properties related to the 'electron', an elementary particle whose 
radius is of the order of 10~ 16 cm. Indeed, as the spin function is arbitrary, there is no loss in 
generality, even if we assume it to be almost a constant (that is, the quantized value of the spin 
of the electron). 

6. THE NEGATIVE ENERGY DENSITY MODEL 

Let us have a closer observation of the results of the previous section 4. The equation (43) related 
to matter-energy density has the spin density part in the first term where spin density is defined 
as s = 35/47ra 3 , where S is the spin of electron the quantized value of which is S = h/2. Then 
substituting the standard values for different parameters in the relativistic units, as mentioned 
in the introductory part, the numerical value for the matter-energy density (43) can be shown 
as 

p = 6.14 x 1(T 37 - 6.81 x 10 27 (1(T 32 - r 2 ). (51) 

The first term related to spin, being of the order of 10~ 37 , is too small compared to the last term 
and hence the spin contribution is negligible. Now, the equation (51) indicates that the central 
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density at r = is negative and its magnitude is about 10 -5 . On the other hand, the total 
density at the boundary, r = a, is positive as usual with the numerical value about 10 -37 . This 
change in the sign of the energy density is because of the presence of the spin term in equation 
(52) which, indeed, is the contribution of the EC theory. In the absence of spin, however, we 
could have negative and zero densities at the centre and boundary of the electron respectively. 
This change in the sign again indicates that the central negative value gradually increases along 
the radius and somewhere, in the region < r c < a, it becomes zero, where r c is the critical 
radius. Obviously, the amount of negative energy density is less than its positive counterpart 
the balance of which ultimately provides the net density as the positive one. 
It is already mentioned that though, in general, for any spherical fluid distribution the density 
on the surface should be zero we are getting here some non-zero value for it. This finite value 
is solely coming from the spin contributed part 2tts 2 . Thus for the vanishing spin the situation 
corresponds to the general behaviour (vide equation (17) of Herrera and Varela, 1994). In this 
context it is also to be noted here that up to the critical radius behaviour of our model is similar 
to that of Herrera and Varela [1994]. Beyond this cut off radius the energy density is regulated 
by spin which makes the overall density of the model as positive. This particular aspect lack in 
the model of Herrera and Varela [1994] where the total energy density is a negative quantity. 
This increase of matter-energy density due to spin density can probably be accounted for the 
kinetic energy through the angular motion of the electron here. 

Similar kind of examination is also possible for the pressures, radial and tangential, both. The 
radial pressure, in this case, takes the form as 

p = 6.14 x 1(T 37 + 6.81 x 10 27 (1(T 32 - r 2 ). (52) 

However, pressure is throughout positive here from the centre to boundary. These results, i.e. 
negative energy density and positive pressure, are in accordance with the pure charge condition 
(37) which reads as p = — p r + Atts 2 . 

7. CONCLUSIONS 

(i) The possible origin of the intriguing negative matter-energy density in the work of Cooper- 
stock and Rosen [1989], Bonnor and Cooperstock [1989], Herrera and Varela [1994] and present 
paper may be due to the finiteness of the total mass of the Reissner-Nordstrom solution [Papa- 
petrou, 1974; Visser, 1989]. Since the electrostatic energy of a point charge is infinite, the only 
way to produce a finite total mass is the presence of an infinite amount of negative energy at the 



9 



center of symmetry. According to Bonnor and Cooperstock [1989] the negativity of the energy 
density and hence the active gravitational mass is consistent with the phenomenon, known as the 
Reissner-Nordstrom repulsion [de la Cruz and Israel, 1967; Cohen and Gautreau, 1979; Tiwari 
et al., 1984; Cooperstock and Rosen, 1989]. In this regards Bonnor and Cooperstock [1989] also 
have discussed about the singularity theorems of general relativity [Hawking and Ellis, 1973]. 
They have shown that the negative regions are liable to exit over distances of order 10~ 13 cm 
and as the proof of the singularity theorems depends on the manifold structure of space-time 
valid down to lengths of order 10 -15 cm so might break down below this. On the other hand, in 
the context of Einstein-Cartan theory of gravitation the idea of negative mass is not a new one 
as stated by Sabbata and Sivaram: "...torsion provides a natural framework for the description 
the negative mass under extreme conditions of such as in the early universe, when a transition 
from positive to negative mass can take place." 

(ii) We have considered in the present paper an extended static spherically symmetric distribu- 
tion of an elementary particle like electron having the radius of the order of 1CT 16 cm, and even 
if for the finite size of the physical system the spin in Einstein-Cartan theory can be related to 
orbital rotation (which indeed is not the case) , for systems of dimensions of subatomic particle 
the orbital rotation loses its meaning. In this case, the only way is to take the spin to be the 
'intrinsic angular momentum', that is, the spin of quantum mechanical origin (in our problem 
since s is arbitrary, we can consider its quantized value or an average value). In this respect 
we would like to quote here from Hehl et al. [1974], "It is crucial to note that spin in U4 theory 
is canonical spin, that is, the intrinsic spin of elementary particles, not the so-called spin of 
galaxies or planets." 

(iii) Following other authors [Prasanna, 1975; Raychaudhuri, 1979; Tiwari and Ray, 1997], in 
the present work the spins of all the individual particles are assumed to be oriented along the 
radial axis of the spherical systems. As to how this alignment is brought about is not very much 
clear. We have discussed here only a few possible ways of realizing this situation. According to 
Raychaudhuri [1979], in general, there will be an interaction between the spins of the particles 
and the magnetic field. The overall effect is the alignment of the spins. In this context, Prasanna 
[1975] mentions that such an alignment may be meaningful either in the case of spherical sym- 
metry when magnetic field is present or else one has to consider axially symmetric field. As 
stated above, our view point is that in the case of physical systems of the size of the electron, 
the radial alignment of the spin is not ruled out. The solution obtained supports this view. 
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(iv) Though our present approach via Einstein-Cartan theory to inject spin may be interesting, 
we feel even that there should have some room to discuss the relationship of our work with an 
alternative means to provide spin and magnetic moment. This is, we think, possible through 
Dirac-Maxwell theory where spin and magnetic moment are naturally incorporated through the 
Dirac spin [Bohun and Cooperstock, 1999; Lisi, 1995] and would like to persue the problem in 
future investigations. 
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